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Abstract 

We show how the two-matrix model and Toda lattice hierarchy presented in a previous 
paper can be solved exactly: we obtain compact formulas for correlators of pure tachyonic 
states at every genus. We then extend the model to incorporate a set of discrete states 
organized in finite dimensional s/2 representations. We solve also this extended model and 
find the correlators of the discrete states by means of the W constraints and the flow equations. 
Our results coincide with the ones existing in the literature in those cases in which particular 
correlators have been explicitly calculated. We conclude that the extented two-matrix model 
is a realization of the discrete states of c = 1 string theory. 
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1 Introduction 

Matrix models have been in the last few years the object of intense research. This is par- 
ticularly true for constant one-matrix models and time-dependent one-matrix models. Con- 
siderably less attention has been paid to two- and multi-matrix models. We think instead 
two-matrix models have a very rich structure and in fact provide a very convenient descrip- 
tion of what is known as 2D gravity coupled to c < 1 matter. In particular, in a recent paper 
lH] we showed that a two-matrix model underlies the so-called c = 1 string theory j^. We 
showed that the correlation functions of the pure discrete tachyonic states of the latter can 
be found in the former and that the underlying integrable hierarchy is the same for both. In 
|||] we did not discuss however the full spectrum of discrete states ||3|,0|||5|. As is well-known, 
in the c = 1 string theory there appear discrete states labeled, say, by integers r and s (see 
below) and organized according to sl2 finite dimensional representations; the cases r = or 
s = correspond to the pure discrete tachyons mentioned above. Occasionally we will call 
the discrete states with r 7^ and s 7^ simply extra states. 

The present paper is a completion of Q. On the one hand we show that the two-matrix 
model of that reference constitutes a powerful tool to compute correlation functions in higher 
genera: we exhibit in particular compact formulas for correlators at arbitrary genus; these 
formulas are calculated in an amazingly simple way. On the other hand we show that the 
two-matrix model can accommodate a full set of discrete states Xr,s organized according to 
the finite dimensional representations of s/2, exactly like the discrete states of the c = 1 string 
theory. To this end we have to modify the two-matrix model discussed in |]l|] by adding new 
interaction terms (like in |^]). The VK-constraints of the modified model turn out to be a 
very efficient tool to calculate the correlation functions of the discrete states at any genus. 
Whenever corresponding results in c = 1 string theory are available in the literature we show 
that they coincide with ours - but our results are much more general. As a consequence we 
consider our modified two-matrix model as a realization of the discrete states of c = 1 string 
theory. 

The paper is organized as follows. In section 2 we present once again the model of [Q] and 
define its extension. We then present the corresponding integrable hierarchy (the extended 
Toda lattice hierarchy) and the corresponding W constraints. In section 3 we define our dis- 
crete states and their s/2 representation properties. We then start the program of calculating 
their correlation functions (CF's) in a definite small phase space, defined at the beginning of 
section 3.2: after some general theorems we concentrate on genus correlators and give a few 
compact formulas. Section 4 is devoted to a discussion of a dispersionless Toda hierarchy, 
which underlies the previous calculations. In section 5 we show the connection with the Pen- 
ner model. In section 6 we return to CF's in higher genera and derive the compact formulas 
announced above. In section 7 we discuss the model in a larger phase space. Appendix A 
is devoted to the discussion of the subtleties of the continuum limit, while in Appendix B 
we collect a few elements in favor of a topological field theory interpretation of c = 1 string 
theory. 



2 The extended two— matrix model 

The model of two Hermitean N x N matrices Mi and M2 , considered in [^ , Q , was introduced 
in terms of the partition function 

ZN{t, c) = j dMidMse^"^, U = Vi + V2+ 5M1M2 (2.1) 

with potentials 

00 
K. = EWK a = 1,2. (2.2) 

r=l 

Although nothing prevent us from solving this model with A^ finite, the spirit of this paper 
is to eventually take N ^ 00; the renormalized N will be called x and identified with the 
cosmological constant. 

After integration of the angular variables eq.( |2.lD can be written as 

N 



ZN{t,C) = const /"fj [|dAc,,iA(Ai)A(A2)exp[/, (2.3) 

•^ a=li=l 

00 N 00 N N 

U = E*i-'tE^M + EvE^2,+5E^M^2,i (2-4) 

fc=l j=l r=l j=l j=l 

where A stands for the Vandermonde determinant. We refer to this model as the ordinary 
two-matrix model. This model can only accommodate the pure tachyonic states. If we want 
to describe the other discrete states we have to enlarge the model in the same way as we did 
in section 4 of [P, that is we add new interaction terms as follows 

00 N 00 N N 

U^U = Y, tl,fc E ^1. + E *2,r E ^2, + V, V = E 9k,r E Al.^2, 

fc=l j=l r=l 1=1 k,r>l i=l 

where gi^i = g. We notice that this corresponds, in general, to adding to the original model 
( |2.1|) not terms Tr{MiM2), but rather composite operators of the form Tr{DiD2), where Z?i 
and D2 are the diagonal eigenvalue matrices of Mi and M2, respectively. 

We will refer to this new model as to the extended two-matrix model. We can generalize 
to this model everything we did for the original one ( |2.lD . In particular we can map this 
functional integral problem into a linear integrable system together with definite coupling 
constraints. In the following we will not repeat the procedure (see [0]), but simply write down 
definitions and results. First let us remind our conventions. For any matrix M, we define 

Mi{j) = Mj^j.i. 







iM),,= 


^-4 


Mij = 


M,„ 


As usual 


we 


introduce the natural 


gradation 












deg[Eij] 


= j- 



and, for any given matrix M, if all its non-zero elements have degrees in the interval [a, 6], 
then we will simply write: M E [a, h\. Moreover M_|_ will denote the upper triangular part of 



M (including the main diagonal), while M_ = M - Mj^. We will write Tr(M) = Eilo^ ^a- 
Finally, for the sake of compactness, we will often use the convention 

As usual we introduce the orthogonal polynomials 

^n(-^i) = A" + lower powers, i]n{^2) = ^2 + lower powers 

and the orthogonal functions 

^n(Ai) = e^^(^i)e„(Ai), <I>„(A2) = e^^(^2)7?„(A2). 
The orthogonality relations are 

^ dXidX2^n{Xl)e^^^''^'^^rrX>^2) = Snmhn{t,g). (2.5) 

In these equations 

oo 

K.(Aa) = ^WA^, a = l,2,, V(Ai,A2)= ^ <7fc,rA^^ 

r=l k,r>l 

Using these relations and the properties of the Vandermonde determinants, one can easily 
calculate the partition function 

Af-l 

ZN{t,g) = const A^! J| hi (2.6) 

i=0 

Knowing the partition function means knowing the coefficients hn{t,g). 

We will denote the semi-infinite column vectors with components ^q,'^i,^2, ■ ■ ■ , and 
$0) ^1) ^2, • • • 5 by ^ and $, respectively. 

Next we introduce the following Q- and P-type matrices 

dAidA2^n(Ai)A„e^(^i'^^)$„(A2) = Qnmia)h^, a = 1,2. (2.7a) 

dAidA2(^^n(Al))e^(^i'^2)$™(A2) = Pnm{l)hm (2.7b) 

dAi(iA2^„(Ai)e^(^i'^2)(— $™(A2))= P„,n{2)K (2.7c) 

Both (5(1) and Q(2) are Jacobi matrices: their pure upper triangular part is /+ = J2i Ei,i+i- 
Here come now the three basic elements of our analysis. The first is provided by the 
following coupling constraints 

P(l) + Y. rgrM^Y'^QC^r = 0, Vi2) + ^ sgr,sQilYQi2r-' = (2.8) 

r,s>l '',.5>1 

These coupling conditions lead to the VFi+oo~constraints (see below). From them it follows at 
once that 

g(l)G [-cx,,l], Q(2)G [-l,cx)] 



The second element are the discrete hnear systems associated to the extended two-matrix 
model. The first discrete linear system is 



f Q(l)^(Ai) = Ai^(Ai), 



r > 1 



dgr,. 



■^(Ai 



Q'-(1)Q^(2) M/(Ai 



r > 0, s > 1 



[ |^*(Ai)=P(l)M/(Ai). 



The corresponding consistency conditions are 

d 

dgrfi 
d 



dgr 



P{l) = [Q\{l),P{l)l 
P(l) = [P(l),fQ'-(l)Q'(2 



(2.9) 





(2.10a) 


r > 0, s > 1, r + s > 1 


(2.10b) 


r > 1 


(2.10c) 


r > 0, s > 1 


(2.10d) 



The second discrete linear system is 
f Q(2)$(A2) = A2$(A2), 
9 .$(A2) = Q^(2)$(A2) 



d9o,s 



S>1 



dgr 



:^(A2 



Q«(2)Q-(1) $(A2 



[ |^<I>(A2) = P(2)v&(A2). 



The corresponding consistency conditions are 



Q(2), P(2)]=l 
^ Q(2) = [(q'-(1)Q^(2)) , Q(2)], 
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dgr,s 
d 

dgo,s 



P(2) = [P(2),(Q^(2)Q^(1) 
P(2) = [Q;(2),P(2)], s>l 



r > l,s > 



(2.11) 





(2.12a) 


r>l, s>0, r + s>l 


(2.12b) 


> l,s > 


(2.12c) 




(2.12d) 



The third ingredient we need is the hnk between the quantities that appear in the hnear 
systems and in the couphng conditions with the original partition function. We have 



d 



■lnZN{t,g)=TT[Q-il)Q^{2)], 



(2.13) 



It is evident that, by using the above flow equations, we can express ah the derivatives of Zjv 
in terms of the elements of the Q matrices. For example 

^ lnZ^(t,<7) = (Q^(a))^^^_, a = 1,2 (2.14) 



dti^ldta,r ^ 'JN,N-1 

Knowing all the derivatives with respect to the coupling parameters we can reconstruct the 
partition function up to an overall integration constant (depending only on A^). 

The three elements just introduced provide a larger definition of our system. In fact we 
notice that the initial path integral exists only in a very restricted region of the parameter 
space. The above formulas, whenever they are meaningful, provide the continuation of our 
system to regions of the parameter space where the path-integral is ill-defined. 

We will use the following coordinatization for Q{'i) and Q{2) 

oo oo 

g(l) = /+ + ^^az(i)ii;i,i_i, Q{2) = 1+ + Y^Y.^i{i)Ei^i_i (2.15) 

i 1=0 i «=0 

One can immediately see that, for example, 

Q+(l)). =5j,i+i + ao(i)5ij, (Q^{2)),=Ri5j.i.i (2.16) 



where i?i+i = /ij+i//ij. As a consequence of this coordinatization, eg. ( 2. 14 ) gives in particular 
the important relation 

\xiZN{t,g)=RN (2.17) 



dtlAdt2A 



and 



92 

--— — In Rj = Rj+i - 2R, + Rj^i (2.18) 

dti^idt2,i 

In our model, the four matrices Q{1),Q{2) and P(1),P(2) are not totally independent. In 
fact the two P matrices can be expressed in terms of the two Q Jacobi matrices as follows 

Pnmil) =Y.^h,rQ'n-Ji^) +nQ~i{l) + Y.^^QnL-H^y, (2-19) 

r=l r=l ^'^ 

Pnm{2) = E rt2,.Q;-^(2) + mQ-ii2) + ^ Q--i(2)^-^. (2.20) 

r=l r=l '^ 

It is worth proving these formulas in some detail. Let us see the first equation as an example. 
From the spectral equation in the first discrete linear system, we see that the polynomials 
Cn(Ai) satisfy the recursion relation 

n+l 



AlCn(Al) = 2^ <5(l)nm'?m(Al) 



m=0 

Solving this recursion relation we get 



Cn(Ai) = det(Ai - g(l) 



The determinant on the RHS refers to the principal minor dehmited by (and including) the 
n-th row and n-th column. Equivalently 

ln^„(Ai) = nlnAi -^ 



^ /A*: dti . 

■r=l J- -^' 



This implies 



In ^„(Ai) = ^ ti^A\ + n In Ai - ^ 



1 dlnZn 



, -; r\\ dti r 

r=l r=l ^ ' 

Taking the derivative with respect to Ai, and using the spectral relation, we have 



(2.21) 



dXi 



r=l r=l ' 



^n(Al) 



ainz„ 



(E ^h,rQ'nn!i^) + ^Qn^(l) + E -^QnL'Hl))^m{Xl 



r=l r=l 

From the last equality we read off formula ( p. 19 ). Similarly we obtain 

1 ainZ„ 

z — ^ ^ 



In $„(A2) = E *2,r A2 + n In A2 - E 



— rX^2 dt2,r 



(2.22) 



from which we can extract ( p. 20 ). 

Among the two sets of flow equations eqs.(2.10b-2.10d) and eqs.( p.l2b -2.12d), we can 
choose as independent the following two 



d 



-Q(1) = [Q(1), (g^(l)Q^(2))j, 
-Q{2) = [{q^{1)Q\2)) , , g(2)]. 



(2.23a) 
(2.23b) 



It is easy to see that if we switch off the new extra couplings (i.e. we set 5^,8 to zero except 
50,^5 5r,o S'lid 5i,i), we recover the usual 2-dimensional Toda lattice hierarchy. Therefore the 
above equations are an extension of Toda lattice hierarchy. Their integrability is guaranteed 
by the commutativity of the flows. 



2.1 Wi+00 constraints 

The Wi+00 constraints (or simply 1^-constraints) on the partition function for the extended 
two-matrix model were obtained in |g]. They take the form 



(£M(1) - {-irT^\l))ZN{t-g) =0, r > 1; n>-r. 

(£M(2) - {-lYT^H2))ZN{t-g) =0, r > 1, n>-r 



(2.24a) 
(2.24b) 



[r] 



H/ 



The generators £« (1) are differential operators involving N and ti^^, while £« (2) have 
the same form with ii^ replaced by t2,A;- The T^ (1) are differential operators involving the 



couplings Qr^s (except grfl) and commute with the £« (1). The Tn (2) are differential operators 
involving the couplings gr,s (except (7o,s) and commute with the £« (2). 
The vCn (1) satisfy a H^i+oo algebra 

[C^:\l\d^il)] = {sn-rm)C^:^:^^'\l) + ..., (2.25) 

for r, s > 1; n > —r,m > —s. Here dots denote lower than r + s — 1 rank operators. The 
algebra of the £« (2) is just a copy of the above one, and is isomorphic to the algebra satisfied 
by the Ti[\l) and by the TJ{\2) separately. 

There is a compact way to write down the above generators. It consists in introducing the 
current 



oo AT ^^^ Pi 

J{z) = J2 rtrz'-' + - + E ^"'" V (2.26) 

r=l r=l 

and defining the density 

^^''^(') = ^--{-9 + Ji^)Y^'---^ (2-27) 

Then £^ can be recovered as 

4"l = Res,=o(^'"l(^)^"+'') (2.28) 

The above definition holds for both the 1 and 2 sector. As for the TA operators we proceed 
as follows. We define the bifocal densities 

Hi{z,w) = Y^ rgr^sZ^'^^w"'^, H2{z,w) = ^ sgr^sz'^^w^'^ , 

r,s>0 r,s>0 

K{z,w)= E z-'-'^w-'-^ ^ 



r,s>0 ^^^'* 

where r and s are never simultaneously 0, and 

Wfl =: K(w-^d, + Hiy : -1, wf' =: K(^z-^d^ + H2Y : -1 (2.29) 

Then 

tM(1) = Res^=oRes^=o(>vf'z''+"u''') (2.30a) 

tM(2) = Res^=oRes^=o(>vf^z^ti;''+") (2.30b) 

For example 

4^1(1) = Y. r9r,sj-^, 4^1(2) = E ^9r,sjr^—, n > -I 

r s>l Ogr+n,s r s>l '-'Qr.s+n 



3 Discrete states and their correlation functions 

The purpose of this section is to show that the extended two-matrix model can accomodate 
discrete states organized in finite dimensional SI2 representations, analogous to the c = 1 string 
theory ones, and to calculate their correlation functions via the VF-constraints introduced 
before. 

3.1 Definition and properties of the discrete states 

We call discrete states the operators Xr,s coupled to the gr,s- In this definition r = s = 
is excluded. However, for later use we introduce also xo,o = Q as the operator coupled to 
90,0 = ^- From the very definition it is apparent that, classically, Xr,s is represented by 
J2i=i Kl^ti where, to simplify our notation we set Ai = A and A2 = /i. As a consequence these 
states carry a built-in SI2 structure. For let us define 

lA/. d d ^ ^ J^. d ^ ^ d 



Then 



^=2g('^5A--^^a^)' ""-^^g'^a^II' ''-^S'^^^ 



[H,E±] = ±E±, [E+,E^] = 2H 



and 



HXr,s = -^{r - s)Xr,s, E+Xr,s = SXr+l,s-l, E_Xr,s = rXr-l,s+l 

Therefore the set {Xr,s = Si=i Kl^t^ r + s = n} form an (unnormalized) representation of 
this algebra of dimension n + 1. 

We can do even better and introduce the new states 



{r + s] 



^r,s = \i —nr^' 



risi 



r,s 



endowed with the scalar product 

One can easily verify that 

^r,s = \j,^TT' >, r = j + m, s = j — m 

is the standard basis for finite dimensional SI2 representations. It is also easy to see that the 
products Xr,sXr',s' with r + s = 2j, r' + s' = 2j' span the tensor product of the representations 
j and j' . It can be therefore decomposed into irreducible representations 



^ri,si^r2,S2 



^C{ri,si,r2,S2\r,s)uJr,s (3.1) 



where the summation extends to the range M(ri — S2,r2 — si) < r < ri + r2 and M(si — 
r2,S2 — ri) < s < si + S2- M{x,y) means the maximum between x and y. Here C are the 

9 



standard Clebsh-Gordan coefficients expressed in terms of the labels r, s instead of j, m. A 
specification is in order: when j = j' the symmetric representations of the Clebsh-Gordan 
series are absent in these classical products. 

To end this introduction let us recall that Xrfi and xo,s coincide with the operators t^ 
and (jg introduced in ||l|, respectively, which, in turn, were identified with the purely tachyonic 
states % and T_s of the c = 1 string theory. It is therefore natural to try to interpret the Xr,s 
as representatives of the discrete states of the c = 1 string theory. 

3.2 Correlation functions of the discrete states. Generalities 

The correlation functions of the extended two-matrix model are in general defined by 

Our main purpose in this paper is to calculate the correlation functions in a very simple 
small phase space Sq: we set gij = except for gi^i = g, which is left undetermined (but in 
the topological application we will set g = —1). As a consequence the CF's will be functions 
of g and A^. We denote by the symbol < • > the CF's calculated is such a small phase space 

0- 

In order to calculate the correlation functions we have simply to write down the W con- 
straints in terms of them. We obtain a set of (overdetermined) algebraic equations which in 
general one can solve recursively. Let us see first of all two very general lemmas. 

Lemma 1. The CF satisfy the symmetry property 

< Xn,si ■ ■ ■ Xr^,s„ > = < Xsuri ■ ■ ■ Xs„,r„ > (3.2) 



This is due to the symmetry of the W constraints ( |2.24a| ) and ( ^.24b ) and of the small phase 
space under the exchange 1 ^^ 2. 

Lemma 2. In the above defined small phase space the CF's satisfy 

{ri + . . . + rn - si - . . . - Sn) < Xn,si ■ ■ ■ Xr^,s^ >= (3.3) 



In order to prove this lemma one rewrites the W constraints ( p.24a| ) and (2.24b) with r = 1 
and n = as follows 

Yl ''9r,s < Xr,s > +^N{N + 1) = 
r>l,s>0 

J2 ^9r,s < Xr,s > +1^N{N + 1) = 
r>0,s>l 

respectively. One differentiates these equations w.r.t. gri,si, ■ ■ ■ , grn,s„ and sets gr,s = except 
51,1 = a- One gets 

(ri + . . . + r„) < Xrt,si ■ ■ ■ Xr„,s„ > +g < Xl,lXrt,si ■ ■ ■ Xr,„s„ >= 
(Sl + . . . + Sn) < Xri,si ■ ■ ■ Xr„,s„ > +9 < Xl.lXn.si • • • Xr„,s„ >= 



■fin the following we often use the expressions the model Sq or the critical point So as equivalent to the small 
phase space Sq. 

10 



Subtracting the second equation from the first we obtain the result. 

This lemma reflects the sh structure of the discrete states since it means nothing but the 
conservation of the H eigenvalue (conservation of the third angular momentum component). 

Consider the W constraints 

(41(1) - i-irTl[\l))z^{t;9) =0, r > 1 (3.4) 

and evaluate it in Sq. We get 

£M(1)^JV = ^^i^ + 1) . . . (A^ + r)ZN (3.5) 

tI[\i)Zn = g' < Xr,r > Zn (3.6) 

This allows us to calculate < Xr,r > for any genus. Remember that the genus h and genus 
h — 1 contributions differ by N'^. Therefore the genus h contribution, < Xr,r >h, is given by 

[ r < 2n — 1 

where 

bk{n)= ^ rir2...rfc, hQ{n) = I (3.8) 

l<r\<r2<- ■ ■<r k<n 

Other exact results will be obtained in section 5. From now on in this section we concen- 
trate on calculating genus zero CF's. 

3.3 Correlation functions. Genus results. 

In genus (alias dispersionless limit) one can produce very general compact formulas for the 
CF's. In order to find the W constraints appropriate for genus we proceed as follows (see 
Appendix). We assign a homogeneity degree [•] to each of the involved quantities 

b.,s] = l, [N] = l, [lnZ(0)] = 2 

where the superscript ^^' denotes the genus contribution. Then we keep only the leading 
terms in this degree in the W constraints. Henceforth we denote A'^ by x (for a motivation see 
Appendix A and end of section 4); the latter will be associated later on to the cosmological 
constant. We will denote the genus zero part of every CF < • > by < • >o. 

The simplest result and the first ingredient we need is < Xn,ri >o- This CF we have 



< Xn,n >o= ^ .^ \ (3.9) 



calculated exactly above, (3/7). The genus contribution is 

^n + 1 

For the remaining CF's the relevant W constraints are 

&i{l)Z{t; g, x) = (-I)'TM (l)Z(t; g, x), r > 1 (3.10a) 

Cl\l)Z{t-g,x) = (-l)'rM(l)Z(t;g,x), r > 1, n > 1 (3.10b) 

11 



and the analogous ones with 1^2. But we wiU not have to consider the latter due to Lemma 
1. Performing the degree analysis one easily extracts the dispersionless limit. In particular 
we have 



T}[^ (1) = Yl ^1 • • • ^rgh,n ■ ■ ■ 9ir,jr 



d 



i "1^ , . . . , i J- > 1 

ji,...jV>l 



dOh 



+...+V+n,ii + ...+ir 



(3.11) 



These generators satisfy the algebra 

[rM(l),rW(l)]=(sn-rm)ri';+^-^](l), r,s>l; n > -r, m>-s (3.12) 

which characterizes the area-preserving diffeomorphisms. The algebra of the Tn (2) is just a 
copy of the above one. Similar simplified (but not quite as simple) expressions can be gotten 
for the >C-type generators. 

Now, from ( 3.10aD and ( |3.10b ), respectively, we obtain 



LHS 



22 kl. . . krgki,li ■ ■ ■ gkr,lr < Xkl+... + kr-r,h+-+lr >0 



(3.13) 



ki,....kr>l 
ll,...,lr>l 



and 






I 1 2;'^ 2^ < Xkuo >o • • • < Xfeno >o= 

/ ki+...+ki=n 

i~^y X! ^1 • ■ ■ ^■'■gki,li ■ ■ ■ gkr,lr < Xki+...+kr+n,h + ...+lr >0 

fej,...,fcr>l 
h,...,lr>l 



(3.14) 



The LHS in ( p.l3| ) needs not be written down explicitly (see below). 

The first step consists in differentiating ( 3.13 ) with respect to gn,n-r with n > r and 
evaluating the result in the small phase space Sq. One gets 

a"" < XO,rXn,n-r >0 +rng''~^ < Xn-r,n-r >0= 

since the LHS of (p. 13 ) does not depend on gn,n-r with n > r. Inserting ( |3.9D we obtain 



< XO,rXn,n-r >0- 



-5)" 



r =< Xr,OXn-r,n >0 



(3.15) 



The second equality follows from Lemma 1. Now differentiate ( 3.14| ) w.r.t. gn.n+p with p > 0, 
and evaluate the result in Sq. One gets 



X < Xp,OXn,n+p >0 



^r-1 



{-lyi^rng' " < Xn+p+r-l,n+p+r-l >0 
+ 9 < Xp+r,rXn,n+p >0 



(3.16) 



for only the first terms in the LHS of ( 3.14| ) contributes, due to Lemma 2. Finally, using (|3.15D 
and (|3.9|) one gets 



^ Xp+r,rXn,n+p ^0" 



X 



n+p+r 



{p + r){n + p) 



-g)n+P+r n + p + r 
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"^ Xr,p+rXn+p,n >0 



(3.17) 



One can write down this result in a label-independent way as follows 

^ _ x^ Min,si)Mir2,S2) 

< Xri,siXr2,S2 >0— /_ s.^ y {6. IS) 

where S = ri + r2 = si + S2 and M(r, s) = max{r, s). This formula also holds when the two 
labels of x coincide. 

The procedure just outlined for two-point functions works in general. If we denote by x' 
a generic product of I discrete states, one first obtains < Xrflx'' >o by suitably differentiating 
( |3.13| ), then one derives < Xr+p,rx'' >o by suitably differentiating ( 3. 141 ). For the generic three 
point functions one gets 

< Xrj,sjXr2,s2Xr3,s3 >0= j^-^M{ri, Si)M{r2, S2)M{r3, S3) (3.19) 

where S = ri + r2 + rs = si + S2 + S3. 

For the n-point functions with n > 3, there is more than one possibility. To explain this 
point we will consider for a while only n-point functions < Xri,si ■ ■ ■ Xr„,s„ >o with r^ / Sk for 
k = 1, . . . ,n (the cases of coincident indices of x ^^^ obtained as limiting cases), and define 
their signature to be {p,q), where p is the number of label r^ > Sk and q is the number of 
labels rfc < Sk- For n-point functions with signature (l,n — 1) we have the general formula 

M{ri,si) . . . M{rn, Sn)(S - 1) . . . (S - n + 3) 

ji + -..+jn (3.20) 

if S > n — 2, and vanishes otherwise. Here we have expressed S in terms of the more standard 
SI2 j labels. 

For other signatures, although there are in principle no obstacles, explicit formulas are 
more elaborate to work out. For example, the 4-point function with signature (2, 2) is given 

by 

< Xrusi ■ ■ ■ Xr4,S4 >0= :E^(n, Si)... M{r4, S4)(S - // - 1) (3.21) 

(-9) 
S = ri + . . . + r4 = Si + . . . + S4, /x = min{\ri - si|, . . . , |r4 - S4I) 

if S > 2, and vanishes otherwise. All the above formulas extend to the boundary where some 
'''k = Sk, and coincide for overlapping values of the labels. 

We remark that the above formulas have been derived for states Xr,s with r and s not 
simultaneously vanishing. To obtain CF's involving p insertions of Q = Xo,0) one has simply 
to differentiate p times with respect to x the corresponding CF without Q insertions. For 
CF's containing only Q insertions, see section 5. 

To conclude this section let us remark that CF's of pure tachyonic states in genus 0, 





^S-n+2 


<Xri,.i.. 


T ^ ^ / 


•A,r„,s„ -U- ^_^^s " 


S = ri + . 


■ ■+rn = Si + ... + S 



calculated with various methods, are numerous in the literature, [^, |^], |1C], pl[| , |1^ , |ll3| , 



[14|. Our results coincide with these up to overall numerical factors. Also formula ( |3.2[1| ) has 
a counterpart in the literature, [p^ ]. 

We can conclude that our results either coincide or extend previous ones about discrete 
states. From this point of view, therefore, our conjecture that Xr,s cLfe representatives of the 
discrete states of c = 1 string theory is confirmed. 
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4 The dispersionless extended Toda lattice hierar- 
chy 

The extended two-matrix model, as we showed in the section 2, is governed by the extended 
Toda lattice hierarchy. In genus zero, the characteristic integrable structure is the disper- 
sionless extended Toda lattice hierarchy. In this section we will give the explicit form of this 
integrable hierarchy and we will show how to obtain, by means of the flow equations and the 
coupling conditions, the same CF's we obtained in the previous section by means of the W 
constraints. The basic ingredient is the continuum counterpart of the relation ( |2.13| ), which 
provides a concrete expression for one-point function, then after solving the coupling con- 
ditions, and making use of the dispersionless Toda lattice hierarchy, we can easily calculate 
all the correlators (besides the CF's among only Q), and recover the results in the previous 
section. 

We proceed as in Q. First we define the continuum quantities in the following way (for 
more details, see Appendix) 



^ ren _ 9r,s 



^ at' ^r.s 



They are the renormalized cosmological constant and couphng constants. Further we define 

Fo(x) = hm ^, C = lim /+ 

where Fq is the genus zero free energy. The second equation is merely symbolic and simply 
means that C is the continuum counterpart of /+. If we define a matrix p = J2n ''^^n,n, it is 
easy to see that we have 

[I+,P]=U, (4.1) 



How this equation enters the Toda lattice hierarchy is explained, for example, in ||15[]. The 
continuum counterpart gives the following basic Poisson bracket 

{C,x} = C. (4.2) 

This allows us to establish the following correspondence 

AT [ , ] ^ { , } (4.3) 



and, similarly, 

^Tr = ^ E ^ L dx (4.4) 



1 „ 1 ^-^ 



N N „ 







together with the replacements 



Q(l) ^ L, Q{2) ^ L^ 
P(l) -^ M, P(2) -^ M 
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where 



R 



L = c + E«^c-', £ = y + E^c 



1=0 



c 



1=0 



Ri 



(4.5) 



Here a; and bi are the continuum fields that replace the lattice fields ai{i) and bi{i) of eq.( 2.15 ). 

We stress that the above replacements holds only in genus 0. The relation (4.3) and (4.4) 
suggest that matrix size N plays the role of Plank constant in a Euclidean theory. In fact 
we need a Wick rotation N — > l/{ih), in order to be able to compare our results with those 
obtained from a Minkowski space theory. 

Next we define the operation |j 



and eqs.( 2.1£ , 2.20 ) become 

M -- 



R 



^if) = f, V function / 



.=1 r=l ^^1.- 



a(M) = £ rt2,rL'-' + xL'' + £ ^Z"^"! 



r=l 



rtl 5*2, 



(4.6) 

(4.7) 
(4.8) 



After this preparation, we can make the following statements: 
The dispersionless version of the coupling conditions (|^.4) is 



(4.9) 



M+ Y. rgr,sL^'^L' = 0, a{M) + ^ sgr,sL''L'-^ = 0, 

T/ie dispersionless limit of the extended 2-dimensional Toda lattice integrable hierarchy ( 2.10b 

dL 



2.12b) is 



ren 

r,s 



dg 
dL 



{L, {nn-}, 
{{L^Ln+, L}. 



(4.10a) 
(4.10b) 



Here the subscript + denotes the part containing non-negative powers of C, while ~ indicates 
the complementary part. The commutativity of the flows directly follows from its discrete 
version. Parallel hierarchies can be obtained by applying the a operator to these equations. 

Next, the continuum version of ( 2.13| ) provides the link between the free energy and Lax 
operators, i.e. 



d 



nven 
■Jr,s 



-F 



{UL'\,){y)dy 



(4.11) 



where subscript '(q)' means that we select the coefficient of zero-th power term of ("• We 
emphasize that this formula is valid on the full phase space. It opens a way to calculate CF's 



^Notice that, with respect to 
their role. 



we have slightly changed our notation, precisely L and d{L) have interchanged 
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by simply differentiating both sides with respect to the appropriate couphng constants and 
use eqs.(4.10a, 4.10b| ). Therefore this equation together with the integrable hierarchy and the 



coupling conditions completely determines the correlators in genus zero. 

Let us see this formalism at work. As a simple example we compute the CF's in the small 
phase space Sq. In such a case, from eqs.(4.9, 4.7,4^), we find 



C. L^^ 



Therefore 








dF 

< X.. >0= ^^ren = 


/^(L^L^)(0)(.V.^...;^^(4 


Then, from 








dg'rT 


(sx'-\'-'+H-g)-', s>r 
"^ ^ 1 0, s<r 



and 



(4.12) 



^^ -{{L^Zn^, l^^l^-'C-'^ya)-^-'' -^- (4.13) 



dg'^"^ Lv /-r> J j^ Q^ ^ ^ ^ 



we obtain 



d'^F _M{r,s)M{k,l) x'^^^ 

< Xr,sXk,l >0- QgrcnQgren " ^7^^^^ ^_^y_,j,dr+k,s+l 



Q^p ^r+k+n—1 

Jr,s '^9kl ^9ri^ 



< Xr,sXk,lXn,m >0= QgrenQgrcnQgren^ = M{r, s)M{k,l)M{n,m)j—^^:^j^Or+k+n,s+l+m 



and so on. I.e. we recover the results of the previous section []. 

This is perhaps the right place to make a comment on the relation between the method 
used in this section and the one used in the previous section. In the previous section we 
used the W constraints to calculate CF's. Generally speaking the W constraints contain all 
the available information of a system, if we exclude the N- dependent constant mentioned 
after eq.( 2.14| ). The W constraints were derived by using both flow equations and coupling 



conditions. In this section we have been using both flow equations and coupling conditions. 
So it is not surprising that we find the same results. The only difference is that in this section 
we have been using a continuum formalism, while in the previous one we used the discrete 
formalism. In other words, in the present section we have first taken N to infinity (with the 
necessary rescalings and renormalizations) and carried out the calculations afterwards. In the 
previous section we first carried out the calculations with N finite: taking A^ to infinity in the 
final results (with the necessary rescalings) leads to the same renormalized results as in this 
section - replacing A^ with the symbol x as we did in the last section understands the passage 
from discrete quantities to (renormalized) continuum ones. 



^The suggestion to represent the discrete states by means of monomials of L and L is also present in ref.llf 
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5 Connection with the Penner model 

The model we studied in the previous sections, characterized by the small phase space Sq, 
has a topological nature. Remarkable evidence of this is provided by its connection with the 
Penner model, ||2l[|. Such connection is to be found in the famous A^-dependent integration 
constant (not determined by the W constraints) which was mentioned after eq.( ^.14 ). 
From eg. (p .171) and eq.(|2.18|) we obtain 



^2 ^2 

In Rn = j—^r—{Fn+l - 2F„ + F„_i), (5.1) 



Where -F„ = In Zn. After integrating twice with respect to ti^i and ^2,1, we have 

In Rn = Fn+l - 2Fn + Fn-l (5.2) 

We remark that this equation can be directly derived from the expression of the partition 
function (2.6). Possible integration constants have been dropped since they are irrelevant for 



the present problem (they provide analytic contributions to the free energy). 

At this point it is more convenient to shift to a continuum formalism (see Appendix) by 
introducing the continuum quantities 

x^^, 6^-, F(x,e) = ^, (5.3) 

as we did in the previous section (from now on the quantities we mention are understood to 
be the renormalized ones). However now we are interested not only in genus 0, but also in 
the higher genera. As is well-known the contributions to the free energy corresponding to 
different genera rescale in different ways. The appropriate genus expansion for the free energy 
is as follows 

oo 

F{x,e) = Y,Fu(x)^'^ (5.4) 

where F^ is the h genus contribution. We also set 

Thus the continuum version of the equation (|5.2| ) is 

\nR(x) = -^ eyiY>{edx) +exp(— eS^) — 2 F{x,e). (5.5) 

e^ L J 

Once we know R{x), we can solve this equation to obtain the free energy. 

Now at the particular critical point {g = —1) we are considering, R(x) = x, which corre- 
sponds to the genus contribution, while higher genus contributions of R vanish (see below, 
section 6). Therefore we can easily obtain the recursion relations among the free energy in 
different genera 

dlFo{x) = lnR{x) 
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and 

n , 
J2j^€Fn-iix) = 0, Vn>l. (5.6) 

Using these recursion relations we can obtain the free energy for any genus. For example, 

Fiix) = -j^^^x, ^2(2;) = -—x^^ 

and so on. In general we have 

F,^(x) = —^^^x'-^\ V/i > 2. (5.7) 

^ ^ 2h{2h-2) ' - ^ ^ 

This assertion can be very easily proved by means of the properties of the Bernoulli numbers 
B21, in particular 

^^ B21 _ 1 

^ 2/(2/-2)!(2n-2/)! (2n)! 

Eq.( |5.7| ) is the free energy of the Penner model. We see that Fh{x){h > 2) indeed have scales 
according to the power 2(1 — h) while Fq and Fi exhibit logarithmic scaling violation, which 
is a typical feature of c = 1 string theory coupled to gravity. The coefficients of the powers of 
X are topological numbers. We recall that the (virtual) Euler characteristics of moduli space 
of Riemann surfaces with n punctures in genus h was computed by Harer and Zagier |2C] and 
rederived by Penner, [pi|], with Feynman graph techniques 

(„)_ {-ir(2h-3 + n)l{2h-l)^ 

Therefore 

Fh = xTx^-^\ < g- >,= nlxPx'-"^-^ (5.8) 

In particular we see that Q, the operator coupled to x, is to be interpreted as the puncture 
operator in the present theory [ p2[ . 

It is a common belief that c = 1 string theory is a topological field theory with primaries, 
puncture equation, recursion relations and a Landau-Ginzburg formulation |16|,|17|,|p^,|19|. 



On the basis of our results, which are far more complete than the ones preexisting in the 
literature, we have reconstructed many sparse elements of such conjecture but no compelling 
evidence. For this reason we do not discuss them here, but collect them in Appendix B. 

6 Correlation functions. Higher genera. 

To calculate correlators in higher genus we can use the VK-constraints as we did for genus 
in section 3. However this method becomes more and more cumbersome as we pass to 
higher and higher genera. Just for the sake of comparison we will see an example below. As 
was anticipated in section 4, however, a skillful use of the coupling conditions and the flow 
equations may produce surprisingly powerful shortcuts. This is what we are going to show 
next: compact formulas for some correlators at all genera. 
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6.1 Exact correlators at arbitrary genus 

In Sq the Jacobi matrices Q(l) and Q{2) take a very simple form 



OO T 

Q{1)=U, Q(2) = ^ — i?;,z_i (6.1) 



which imphes in particular - compare with ( 2.16| ) 



Ri = — (6.2) 

-9 

Eqs.(|6.1|) can be derived by solving the coupling conditions. In Sq, the coupling conditions 



(2.8) are simply 

P(l) + gQ{2) = 0, P(2) + gQ{l) = 0. (6.3) 

Furthermore, in Sq, eq.( ^.19|) shows us that 

P+(1)=0, Pn,n~l{l)=n 



Combining these equations with the first equality in eqs. (|6.3| ), we get the same Q{2) in 
eqs. ( 16.11 ). Similarly we can obtain Q{1)- In the rest of this section, we will show that eqs.( |6.1[ ), 
together with the extended Toda lattice hierarchy, are sufficient to exactly determine the cor- 
relators. 

6.1.1 The 2-point correlator of purely tachyonic states 

The calculation of the correlation functions consists of the following three steps. 

1. Evaluating the correlator on lattice: Using eq.( ^.13 ) and eqs. (|2.10b| , 2.12b ), we immediately 

get 

^' lnZ,(g)=TY([Q™(l),Q"(2)] 



Evidently it is non- vanishing only if m = n, and a simple calculation shows that 

o2 J i+n-i 

;^j-^l„Z,fe)=Iv([Q»(l),Q»(2)l)= S UR, 



1 ' 



J2 /(/ + l)...(Z + n-l) (6.4) 



2. Shifting to the continuum formalism: The renormalized correlation functions are the deriva- 
tives of F{x) with respect to the renormalized coupling constants 

F{x,e) =< r„r_„ >i.en 



dtf:^dtf:?, 

For simplicity, from now on, we drop the label ren wherever it should appear. From the context 
it will be clear whether the quantities in question are unrenormalized (before the continuum 
limit) or renormalized (after the continuum limit). 
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A straightforward calculation shows 



1 ° 



<TnT^n> = 7 V (2; + /e)(x + (Z + l)e)...(x + (/ + n-l)e) (6.5) 

-, rt— 1 

= T-T^T.(^ + le)ix + {l-l)e)...ix + il-n + l)e) (6.6) 

'^ ^) 1=0 

Separating the contributions from different genera: In order to single out the contributions 
from different genera, we recall once again the topological expansion for the free energy 

oo 

F{x,e) = J2Fhix)e^^ (6.7) 

h=0 

Further we observe that eq.( |6.6[ ) correlator can be re-expressed as 

1 ° 

<r„T_„> = — _-( ^ e-'^^-)x(x + le)...(x + (n-l)e) 



i-g) 



-. 71—1 

(J2e'''-)xix-le)...ix-{n-l)e) (6. 



2h 



= ^ < TnT-n >h e 
h>0 

The second equality exhibits the invariance under change e ^f — e, which amounts to the 
vanishing of the odd powers of e in the last expression. In other words, the CF's have the 
same kind of expansion as the free energy. 
Now let us define two series of numbers 

n-l 

ai{n) = Y,l\ i > 0; 

1=0 

i 

l<h<l2<...<li<n-l p=l 

The first few ai{n) are 

ao{n) = n, ai(n) = ( ^ j , a2(ra) = — - — ( 2 ) ' 

,, n2(n-l)2 ^^ (2n-l)(3n2-3n-l) / n \ 
«3(n) = , a4(n) = (^ 2 j 

The first few /?i(n) are 

Poin) = 1, /3i(n) = ( 2 ) ' ^2(^) = ^^ ( 3 ) ' 

^ , , f n\ f n\ ^,, 15n3 - 30^2 + 5n + 2 / n \ 
^^^") = (^ 2 J (^ 4 J ' ^^(") = 48 [ 5 j 
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Furthermore, from the definition, we directly obtain 

(3n-i{n) = {n- 1)!, A(n) = 0, Vi > n. 

Notice that bi{n) = /3j(n + 1), where bi{n) were defined in subsection 3.2. In terms of these 
numbers, we have the expansion 

n 

x{x + e){x + 2e) . . . {x + (n - l)e) = ^/3i(n)a;""*e* 

i=0 

and 

/=0 1=0 

Plugging these formulas into eq.(|6.8|), and remembering that terms with the odd powers of 
e vanish, we obtain the general two-point function (which include the contributions from all 
the genera) 



1 ^^ ( 

-^ E ^^'^E(-i)M 

^> 0<2h<n 1=0 \ 



I ^ — ^- \„ t„\i3 t„\^n-2h 



<r„r_„>=— — Y^ e^'^^(-iy 2;_J«2.-/(n)A(r^K-^'^ (6.9) 

*- ^' 0<2h<n 1=0 \ / 

where g is the renormalized coupling constant (we dropped the superscript for simplicity). 
The /i-genus correlator is 

n—2h 2/i / , \ 

< T„r_„ >h= ^—— E(-l)' ( 2h~- I j «2/.-/(^)A(n), 2/1 < n; (6.10) 

< T„r_„ >h= 0, 2h > n. 
In particular we have 

< r„r_„ >o= -, ^nx"", (6.11a) 

< Tr^T^n >i= j^n f "" 4 M x^-\ (6.11b) 
^ rj. rj. ^ 1 (3n^ - n - 6)n /^ re + 1 ^ „_4 /^ 1 1 a 

< InJ^-n >2= -, ^ ^ \ X . O.llC 



6.1.2 The 3-point correlators of purely tachyonic states 

Now we turn our attention to the three-point correlation functions. Once again, from the 
Toda lattice hierarchy ( p. 10b , 2.12b| ) and the relation (2.13), we get the correlator on the 
lattice 



J 



Tr([[Q"(l), Qi (2)]+, Q"(l)]) +TY([Q^(1), [Q-(1), QL (2)]]) (6.12) 



dtl^mdti^ndt2,l 
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Therefore in Sq we have to evaluate 

J — 1 i+l i+n i+m i 

Tv([i^, [I™, (/?/_)']]) = si,r,+^Y.{i[ ^p- n ^p- n ^p+ n ^p 

i=0 p=i+l p=i—m+l p=j— n+1 p=i—l+l 

= Si^n+m[ Yl ~ J2 ) II ^P 
i=j—m+l i=j~l+l p=i 



After passing to the continuum, eq.( |6.12 ) becomes 



^l,nA 



m-1 l-l 



< TnT„,T^i >= T^(E e"^'''^ - E e-P''^^)x{x + e){x + 2e)...{x + il- l)e) 

V 9) ^ p=o p=n 

I- m—1 l—l 

- ^''"+I!l(^ eP'^- - E e^^^-)x(x -e)ix-2e)...ix-{l- l)e) 



(-5)'^\t^o 



p=n 



Again, the second equahty confirms that terms with the odd powers of e vanish. Expanding the 
two factors in terms of a, /3 numbers, we get the compact formula for three point correlation 
functions 

l-2h-l 2h / / _ ,■ \ 



i=0 



■(^a2h-i+i{n) + a2h-i+i{m) - a2h-i+i{l)), 2h < I - I; 

< TnTmT^i >h= 0, 2h>l-l. 

In particular for genus one we have 

< TnTmT_i >i= di^n+mj^-y: ^ 1 3 I ^ ^ ^ 

Let us remark, after these two subsections, that whenever two- or three-point functions 
of tachyonic states in genus 1 are available in the literature in explicit form, l^,l^,f^, they 
coincide with the above results up to overall normalization constants. 

6.1.3 Other discrete states 

The previous method does not work only for purely tachyonic states. For example 
dFN 



Tr(Q^(l)Q^(2)) = 5r+sfi{Q''{l)Q'{2) 



d9r,s 

On the other hand, in Sq, 

N-l , N-l 



(q"(1)Q"(2)) = E ^p+1 • • • Rp+n = j^ J^{p+l)...(j) + n) (6.15) 
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The summation can be performed directly, and we obtain 

from which we recover the results of subsection 3.2. Now we can take the continuum limit, 
which amounts to replacing A^ by x, we finally have 

^n-2h+l 

<Xn,n>h=p2h{n+l)^—-— -, 2h<n-l. (6.16) 

(n + l)(-g)" 

In this particular case it is more convenient to carry out the calculation in the discrete and 
take the continuum limit of the results, instead of going to the continuum, as in the previous 
examples, and then carrying out the calculations. The two methods lead to the same results 
as pointed out in the Appendix. 

Let us see another example: the two-point function for x states. 

+ <3(l)'-|(<3(l)"Q(2r)^,(3(2) 

Inserting eq.( |6.lD , we see immediately that this does not vanish only if r + n = s + m. Without 
loss of generality we can suppose n > m and s > r. Then we obtain 

Q2p N N+r-s 

-^ T^ = 2^ Rp+s+m-1 ■ ■ ■ Rp — 2_^ Rp+s-1 ■ ■ ■ RpRp+n-l ■ ■ ■ Rp+n-m (6-17) 

To calculate this correlation function we proceed as in the example just shown. We replace 
the explicit form of Rp, we notice that the genus h contribution in proportional to j\[^+'^-'^h 
and single out the corresponding coefficient. In the relevant formula we take the continuum 
limit, which amounts to replacing N with x, and obtain: 



^ ~X.r,s'X.n,m ^h" 



m+s—2h 
(-5)™+^ 



^^^+^^"^ + ^^ (6.18) 



'^' . (l(s + m-l\. {r-sY+^'--' ( s + m-l + l 

V- B2q_( S + m-l\ ( S + m-l-2q+l\ .2h-l-2q+l 

2^2hm^^^<i \ 2^-1 A 2/.-/-2,+ l }^^-'^ 

for m-\-s> 2h, and zero otherwise. The convention ((— n)!)~^ = is understood when re > 0. 



In deriving the above equation we have made use of a well-known summation formula |23|, 
which involves the Bernoulli numbers -625- Moreover we have introduced one more series of 
numbers 

I 
'yi{s,n,m) = Y^ ^ ^ iii2 . . .ikjk+i ■ ■ -ji (6.19) 

fc=0 0<Ji<i2-..<ife<s— 1 n—m<jf^^i<...<ji<n—l 
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which are non-vanishing only if n > m, and / < tti + s — 1. We have in particular 
7o(n,m;s) = 1 



7i (n,m;s) = —{s{s — l)+m(2n — m—l)) 



72(n, m;s) = —i—2m — 3m +2m +3m +12mn—12m n—\2mn 

9 9 9 9 

+ 1277T, n — 2s + 67Tis + 67n s — 12mns + 9 s 

9 9 9 9 "^ id ' 

— 6ms — 67n s + 12 mns — 10s+3s 



We could of course continue with other more complicated examples. But it is perhaps 
more interesting to make a comparison with the method based on the W constraints. 

6.2 Genus 1 correlators from W constraints 

In this subsection we resume the method of subsection 3.3, to calculate correlation functions 
of discrete states at genus 1. This method is more cumbersome and less efficient than the 
previous one. However it provides an independent check of the previous calculations. We 
proceed as in section 3.3 with an obvious modification: in the expansion according to the 
homogeneity degree (see Appendix), we have to take into account all the contributions up to 
and including the genus we are considering. For example, for genus 1, in the expression for 
Tn (1) we have to take into account the following contributions 

/ d 
T^\'^)= XI ikik2---krgkuii9k2,h---9kr,ira (^-2^) 

fel,fe2,...,fc, C5fci + ...+fe.+n,h+...+/. 

Il^l2,---,lr 
f f) 

+ -kik2 . . . kr-i{ki + ... kr-i - r + l)gki,h ■ ■ ■ 9kr-i,i 



^ Ogkx+...+kr-i+n,h + ...+lr- 

fff 1 ) (J 

H ^ ki . . . kr-2^(,kp - l){kp-2)gk,j, ■ ■ ■ gkr-2,lr-2 



6 p^i - - Ogkj^ + ,,, + kr-2 + n,h + ... + lr-2 

'V { f 1 ) (f 

+ -. kl . . . kr-2 Yl (^P - 1)(^9 - ^)9kuh ■ ■ ■ 9kr^2,lr-2j^ 

4 p,j^i (Jgki+... + kr-2+n,ll + ... + lr-2 

p<q 

The first term in the RHS contains the genus contribution. 

For the sake of conciseness, we will not write down the corresponding expression for £„ (1), 
which is even more complicated. Eq.(6.2C|) gives an idea of the complications we run into by 



proceeding with this method. 

The first step of the calculation consists in evaluating < Xn,n >, which we already know. 
Then by differentiating (C_1{1) — {—1YT_1.jZn = with respect to gn+r,n one obtains 

n+r-2 
< XO,rXn+r,n >1= 24(_ )n+r '^^"' + ^)('^ + ^ ~ l){3n'^ + r^ - n - r - 2) (6.21) 
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Then differentiating f £m(l) — (— l)'"Tm (l))^Ar = with ?ti > 0, with respect to gn,n+m, and 
evaluating the result in So, we obtain 

< Xn,n+mXr+m,r >1= 24(_^)n+m+r ^ + ^ + ^_2 ' ^^'^^^ 



•(—4 + 4771 + 3771 —Am +m +4n + 2mn — 6'm n + 2m n + 7n —12mn 
+ 4m 77 — IOti +677771 +377 + 4r + 27?7r — Gtti r + 27r7 r 

9 9 9 "-t 9 9 

+ 1077r — 247?7nr + IO777 Tir — 1877 r + \Amn r + 6n r + 7r — 127nr 

99 9 9 99 "-t "-i "-tA 

+ 4777 r — 18 77r +lAmnr + 10 ti r — lOr +67nr +677r +3r 
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We have checked that this is exactly equal to eg. (6.1^) above for h = 1, provided of course we 



choose corresponding labels. As a particular case of (|6.22|) or ( |6.21|) , we obtain 

< Xm,OXO,m >1= 24f_ )m ^^ + l)"i^("T' " l)(m - 2) 

which is eq.( |6.11b| ) above. 

In the above calculation actually something new happens with respect to section 3.3: there 
appear 'half-genus' contributions. In Appendix A we explain that F, the free energy, is a sum 
of contributions from different genera, and F^, the genus h contribution to the free energy 
scales like 2(1 — h). This is perfectly satisfactory for the ordinary two-matrix model. But 
for the extended two-matrix we find also contributions that scale like (1 — h), which we call 
with an effective term 'half-genus' contributions. This phenomenon is apparent for example 
in eqs.( |3.5| , p.6| ) and has already been remarked in section 6. Such contributions appear as 
intermediate results in the calculation on CF's via the M^-constraints and are necessary for 
the consistency of the theory. From a practical point of view one has to them into account 
and simply select the integer genus contributions in the final results. As a matter of fact 
such contributions are not new as they appear for example in the intermediate steps of the 
calculation of CF's in models corresponding to n-th KdV hierarchies or other hierarchies. But 
in the present case they appear as contributions to the CF's themselves. They may in fact 
represent a new and interesting phenomenon: they may imply that the extended two-matrix 
model takes contributions not only from compact Riemann surfaces, but also from Riemann 
surfaces with boundary. 

7 Beyond the critical point 

Till now we have only worked in the phase space Sq. What happens if we enlarge it? In this 
section we briefly discuss this point. For the sake of simplicity, we focus our attention on 
genus zero calculations. 

7.1 Perturbing around Sq 

Let us Taylor expand the free energy around Sq. The genus h part takes the form 

Fh{gr,s) = !< h{O0) + 2^ gr,s—K ^ 2^ 9ri,sigr2,S2-K ^ ^••• 

r,s Cl9r,s ri,r2;si,S2 09ri,siOgr2,S2 
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-Fh{So) + ^ gr,r < Xr,r >h (7.1) 

r=l 

/ , "ri+r2,si+S25ri,si5T-2,S2 ^ Xri,siXr2,S2 ^/i + • • • 



^i,»"2;si,s2 



In the second equality, the first term has been given in section 5, the other terms are ob- 
tained by plugging in the correlation functions of the small space Sq. The 'dots' contains 
the contributions from multi-point CF's. The free energy can be divided into two parts: 
coupling-independent and coupling-dependent. These two terms are quite different: the 
coupling-independent term is singular at x = 0, the other is regular; even more crucial the 
coupling-dependent term are completely determined by Wi^oo constraints, while the coupling- 
independent term is not derivable from Wi+oo constraints. In fact, as we have seen in section 
5, it stems from the Toda lattice equation and the coupling conditions. Finally it is just the 
coupling-independent term that give rise to the logarithmic violation of the scaling property 
in genus zero. One should always bear in mind the existence and the difference between these 
two contribution to the partition function. 

7.2 Correlators in S. 

The small phase space S is defined as follows 

S = {gr,s = 0, except gi^i = g, go,r = il,r, 5r,0 = t2,r} 

In other words, we switch off the extra couplings. 

At first we notice that in the continuum limit the genus generators Cn take a very simple 
form. Remember eqs.( p.26 ,2.27, 23^ ) and remark that 



Z-^J{z)Z = J2 rti,rz'-^ + - + Y. 



z 

r=l " r=l 



_ r-l d^O 



which has the same form as M(L). From (3.10b) and ( |3.11| ), we get therefore 



E (n ^rg^.,^-, ^-^ = , },T^^^ iM-^'dL-+\ (7.2) 

n,.S;l>i V=i ^ %i+...+i™+n-mji+...+i„ (n + l)(m + 1) / 

jlv,im>l 

Simply restricting to S we obtain the one-point function 

We remark that this formula is very similar to the period representation of one-point function 



in Landau-Ginzburg representation |24] 



In order to obtain the multi-point functions of the extra states, we differentiate eq.( p'.2[ ) 
with respect to the extra couplings and evaluate the result in S. For example, the two-point 
function of the extra states is 

< Xk,lXn,m >0= . ^.w ^,w ^J^ f M-^+^{L)dL-^+^ < Xn+k-l,ra+l-l >o(7.4) 
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The first term in the RHS means that one should first take the derivative with respect to the 
coupling gk^i (using the dispersionless Toda hierarchy), and then do integration. With the 



terminology of the Liouville approach [25| pH], the first term may be interpreted as "6u//c" 
contribution, while the second term can be viewed as a ^''contacf term of two operators. 
Similarly we have the three-point function 



< Xr,sXk,lXn,m >0 



02 



{n + l){m + l){-g)"'dgr,sdgk 

km 

< Xr,sX 

9 

rkm{m — 1 



- i M"'+\L)dL''+^ 



km rm, 

< Xr,sXn+k~l,m+l-l >0 < Xk,lXn+r-l,m+s~l >0 (7.5) 

9 9 



9' 



< Xn+r+k-2,m+s+l-2 >0 



Without much difficulty, we can derive formulas for the most general CF of extra states on 
the small space 

k k r^ 

^=1 (n + 1) (m + l){-gr ^Ji %m j^ ^ ^ 



^ Y^ Vi)V(2)---V(r )^("^ - 1) . . . (m - r + 1) 

r=l p€Sk 



9' 



(7.6) 



< Xra+ip(r+i)+...+«p(fc)-r',m+jp(r+i) + ---ip(fc)-^' 11 Xip(^s),jp{s) ^0 



s=l 



where p is an element of the symmetric group 5^. The correlation functions involving pure 
tachyonic states are nothing but particular cases of the above formulas. 



7.3 Correlators in Si. 

The above formulas for CF's are implicit. In order to get explicit results we have to solve the 
coupling conditions. This is particularly simple if the small phase space is 



Si = {gr 



0, 



except 51 1 = g, 



h,i, 



t2,l} 



In this case, from the coupling conditions (|4.9| ), we can write down two dispersionless Lax 
operators 



L = C + 



*2,1 



L 



1 ,x . 



together with 



R 



X 



(-9) 



Plugging these expressions into ( [4.8D ,( |7^ and ( |7.4| ), we obtain 



min(n,m) 

< Xn,m >0= 2^ 
r=0 



n 
r 



m, 
r 



r+l j.m~r j.n-r 
^ ''1,1 ''2,1 



(r + l)(-5) 



n+m—r 



(7.7) 
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and 



min(n+fc,m+Z) / i \ / ,7 1 \ r^r j.m+1-r ,n+k-r 

1 X^ n + k-1 \ m + l-1 \x % ^ ^2,1 .„ „x 

< Xk,lXn,m >0=km ^ I ^ _ 1 I I ^ _ 1 I ^(_g)n+m+fc+^-r (^-S) 

/ \/ \/r\/;\ -yj+s .m+l-j-s n+k-i~r 

Y- (■ -^ai- -Ne n\m\k\l\x-' t^^ tg,! 



i /I j j \ r j \ S j ( — Q)n+m+k+l-i-r 

0<i<n,0<j<m \ / \ / \ / \ / \ ^) 

0<r<k,0<s<l 

and so on. As for pure tachyons we find, for example, 



and 



< Tn >o= ^, < T_„ >o= ^ (7.9) 



nm 

< TnTm >0= 0, , < T^nT-m >0= 0, < TnT^m >0= < Xn-l,m-l >0 (7- 10) 

-5 



as well as 



< QTnTm >0— 0, < QT^nT-m >0— 

Tnin{n,m) / -, \ / -, \ ^r -1 j.m-r ,n-r 



r=l 



r-1 /\ r-1 / (_^)n+m-r 



rr rr rr n t^ t^ t^ nml{l — 1) 

< TnTrnTl >o= 0, < TnTmT^l >0= 7 To < Xn+m-2,i~2 >0 

i-ar 

rp rp rp p. rp rp rp ^(^ ~ l)ml 

< T-nT-mT-l >0= 0, < TnT_rnT-l >0= 7 T^ < Xn-2,m+Z-2 >0 • 

These correlators explicitly show that the selection rule implied by Lemma 2, section 3, does 
not exist any more in Si. 

From the above it is evident that the model Si is quite different from Sq. However one 
can still conclude that Q plays the role of puncture operator. Let us see this point in more 
detail. We denote the couplings which are not in the small space Si by gr,s-, and we suppose 
that they are infinitesimal parameters. Then we can evaluate the free energy as follows 

Fh{9r,s) = Fh{Sl) + 2_^gr,s—^ H 2^ 9rusi9r2,S2^ ^ ^■■■ 

r,s '^9r,s ri,r2;si,S2 Cl9ri,siOgr2,S2 

= Fh{So) + ShfihSiX + Y^ gr,s < Xr,s >h+ Yl 5ri,siffr2,S2 < Xri,siXr2,S2 >h + • 

r,s ri,r2;si,S2 

The first two terms are calculated from Toda lattice equation recursively by noting that R = x, 
which shows that Ffi(Si) receives an additional contribution only in genus zero. The remaining 
terms are proportion to the infinitesimal parameters gr,s- Now taking derivatives with respect 
to X, we obtain 

d"'Fh(x) 
< Q" >h= — -K-n \- 6n,iSh,ohsi + g - dependent terms 

Returning to the space Si by setting g^^s = 0, we see that the multi-point functions of 
the operator Q (besides < Q >o ) still represent the Euler characteristics of the punctured 
Riemann surface, i.e. Q still plays the role of puncture operator. 
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8 Conclusion and outlook 

In this paper we have often utilized two different methods to derive the same results. We have 
seen that our results compare very well with the existing ones in the literature. We think 
we have given convincing evidence that the extended two-matrix model not only provides a 
realization of the discrete states of the c = 1 string theory, but is also a very powerful tool for 
solving it. 

As a matter of fact we have not mentioned so far the approach based on the Liouville 
dressing of matter vertex operators, which makes essential use of the Coulomb gas representa- 
tion 1 27 1, 1 28 1, ||2^. Comparisons with this approach and related problems have been discussed 



in the literature we have cited. We do not have anything new to add to them. 

We have not addressed the theme of the ground ring of c = 1 string theory either, |30|. 



We believe a further extension of our model is necessary for us to be able to describe such 
ring. We will return to this elsewhere. 



Appendix A. Discrete and continous formalism. 

In this Appendix we discuss the relation between the discrete and continuous formalism, which 
are both used in this paper, and explain why and when they give coincident results. 

The equations we have to solve to calculate CF's are discrete as far as the lattice variable 
A^ is concerned. We recall that N is the common size of the matrices Mi and M2. As we show 
in the paper one can carry out many calculations sticking to the discrete formalism. This 
is in fact the most natural thing to do. One then takes A — > 00 in the results themselves. 
But shifting to a continuous variable x is often more comfortable and more effective. The 
passage from A to x is regulated by precise rules, which are dictated by the scaling properties 
af the quantities in question. The latter differ genus by genus for the free energy and all the 
quantities connected with it (CF's, fields, etc.). 

In order to split the contributions from different genera we rescale the lattice length. 
Simultaneously we have to rescale all the quantities according to their degree of homogeneity 



[gr,s]=a + pr + js, [N] = a, [Ft] = 2{1 - h)a 

where F^ is the /i-genus contribution of the free energy. Analogous genus expansions hold for 
all the fields ai, bi and R of the theory. The numbers a, f3 and 7 above are arbitrary. For the 
purposes of this paper it is enough to choose a = 1 and /? = 7 = 0. If we denote by n (instead 
of A) the matrix size and imbed the lattice variable n into a finite lattice of size A, then the 
above choice corresponds to rescaling the lattice by A^ so that the lattice spacing becomes 
e = 1/A. We then set x = n/N . This choice has the virtue that the contributions to the free 
energy corresponding to neighboring genera differ by e^. We stick to it throughout the paper. 
After the correct rescalings have been performed, we can consider a; as a continuous variable 
if n and A are large enough and proceed according to the rules of the continuous differential 
calculus. 

Occasionally, however, we avoid such pedantic procedure as to introduce n and A as above, 
and simply say that A^ (the matrix size) is replaced by the renormalized variable x. 
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Two remarks are in order. First, the homogeneity argument is very important in our ap- 



proach, and it is suggested by the Feynman diagram approach of |31|. This means nothing but 
the fact that we look only for solutions with the above homogeneity properties; in particular 
this means that we are interested in genus-expanded solutions. 

Second, passing to the continuous variable x is, in many cases, not strictly necessary for the 
approach to work. As a matter of fact, since in our two-matrix model there is an arbitrarily 
large number of couplings, one should eventually take A^ to infinity: this is particularly evident 
from the calculations of section 6. However one can take the attitude that N is arbitrarily large 
but finite, and one can let A ^ oo (with the necessary rescalings) once all the calculations have 
been carried out. The fact is that the equivalence we establish between fn+i and exp(e9)/(x) 
is (up to a possible rescaling) an exact relation even in the discrete case, at least as long as fn 
is a polynomial in n. As long as the quantities we have to do with are polynomials, the two 
methods (discrete and continuous) give the same results. 

In conclusion, except when transcendental functions are involved, we can follow two 
courses: either we go to the continuum (with the necessary rescalings and renormalizations) 
and then do the calculations, or we carry out the calculations with finite N and take the 
continuum limit (with the appropriate rescalings) on the results. 

It also follows from the above considerations that, although in this paper we have privileged 
the N -^ oo limit, the correlation functions of the two matrix-model studied here can be solved 
with A^ fixed and finite. A model with finite A^ can be of interest in itself. 

Appendix B. Topological field theory properties. 

The so-called c = 1 string theory is believed to be representable as a topological field theory. 
In this Appendix we collect some (non conclusive) evidence of this conjecture. I.e. we try to 
envisage the model Sq studied in section 3 (at g = 1) as a topological field theory model. In 
order to be able to make such assertion one should be able to identify the puncture operator, 
the primary fields and the descendants and show that they define an appropriate metric and 
satisfy a puncture equation and the appropriate recursion relations. 

We have already seen that the puncture operator has to be Q. As primary fields we take 
Tn = Xn,o and r_„ = xo,n, where n is any natural number and Tq = Q, while all the other Xn,m 
are descendants. We recall that r„ and T-n were identified in Q with the purely tachyonic 
states Tn and T_„, respectively, of c = 1 string theory. 

In this Appendix we have to pay some attention to indices. Since, in order to conform 
with the current notation, we have introduced negative labels, we will explicitly point it out 
whenever the labels r, s, ... denote integers and not simply natural numbers. 

The metric (in the topological field theory sense) is given by 

r]k,i =< QTkTi >o (8.1) 

where k and / are integers. The only nonzero elements are (at g = —1) 

d 

Vn~n = V-n^n =< QTnT^n >0= ^ < Xn,OXO,n >0= n^x"~\ r/o,0 = X~^ 

This metric is non-degenerate, the inverse is r] ' with 
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while all the other elements vanish. The associativity condition for the structure constants 
CijM, i,j,k integers, 

is easily seen to be satisfied once we notice that the only nonvanishing three-point functions 
among primaries are 

'^n,m,—n—m — '^—n,—m,n+m — ^ J- —n^ —mJ-n+m ^0 — flTTlyTl -f- TTljX l"--^j 

nm{n — m)x^~^, n>m 
nm(m — n)x"^~^ , n<m 



^n,—m,m—n — '\ / \ m—1 



beside Co^n,m = Vn-n^n+mfl- The primary fields form the commutative associative algebra A 

k I 

where z, j, k, I are integers and Tq = Q. To prove it one has to use 

^ n+m _ "-"^ /-» n _ 1 ^ _ „2 n ^ — 1 

'^n.m — ; ! "-^O.ra — ■■-) '^n.—n — ii' X , (-^O.O — ^ 

n + m 



^ n-m _ I n-m 

where n,m ^0, Cij^ = Cj^i^, and the other structure constants vanish. 

Let us come now to the recursion relations and puncture equations. The recursion relations 
in Sq are 

< Xr,sXrusiXr2,S2 >0= Af (r, s) ^ < Xr-l,s-lTi >o V^''' < TfcXn.siXra.sa >0 (8-3) 

l,k 

where the labels k and I are understood to be integers. The proof is very simple. Suppose for 
example that r > s. Then 

LHS = rM(ri,si)M(r2,S2)x'^+'^i+''2-i 

when r + ri + r2 = s + si + S2 and vanishes otherwise. On the other hand 

RHS = r < Xr^i,s-iTs-r >o 7]'-'^'-' < r,_,Xn,siXr2,s2 >o= rAf(ri,si)M(r2,S2)x''+"i+"^-^ 

when r + ri + r2 = s + si + S2, and vanishes otherwise. The same can be proven for r < s. 
The puncture equations are designed to connect the CF's of of the type 

^ QXr-i,siXr2,S2 ■ ■ ■ Xrn,s„ ^0; 

where the x's are extra states, with CF's including neighboring descendants of them. For 
dimensional reason the latter can only be < Xn, si ••• Xr,-i,sj-i ••• Xr„,s„ >o- A heuristic 
relation which does this job for the CF's (|3.20| ) is the following 

< QXr,,siXr2,S2 ■ ■ ■ Xr„,s„ >0= (8-4) 

^ M{ri,Si) S-1 ^ 

^M(ri-l,s,-l) n <^r^,s,...Xn~i,s^-i...Xr„,s„>o 
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where E = ri + . . . + r„ = si + . . . + s„. In fact the LHS is 

< QXrus^Xr2,s, . ■ ■ Xr„,s„ >o= x^-"+^M(ri, Si) . . . M{rn, s„)(S - 1) . . . (S - n + 2) 
On the other hand the generic term in the RHS of (^^) contains 

< Xri,si ■ ■ ■ Xrfc-l,Sfe-l • • • XTn,Sn >0 = Af (n, Si) . . . M{rk - 1, Sfc - 1) . . . M{rn, s„ 



(S-2)...(S-n + 2)x 



S-n+l 



Summing all the contributions in the RHS of ( |8.4D we obtain the equality with the LHS. 

The trouble with this formula is that it does not extend, as it is, to other CF's and to 
higher genus. Modifications of this formula can be figured out in many cases. For example, 
for one point functions, 

< (1 - e-^)Xr,r >=r < Xr-l,r-l > (8.5) 



is exact and is clearly the generalization of ( |8.4| ) to every genus. But, unfortunately, we do 
not have a general compact formula for all the cases. 

If we accept this point not as an imperfection of our proof but as a peculiarity of the 
puncture operator Q, we can then regard our model as a topological field theory in which 
Q = To is the puncture operator, T^ and T^n {n natural) are the primary fields, while Xn+k,k 
and Xk,n+k, with k positive, are, respectively, the descendants. In particular Xk,k are the 
descendants of Q. 

Before we pass to the LG formulation, a comment about the position g = —1 we made 
throughout this Appendix is in order. The reason for such position is that only for this value 
of the coupling do puncture and recursion equations hold, as is easy to check. Alternatively 
the coupling g can be absorbed into a redefinition of Xr,s'- precisely Xr,s — > {V~9Y~^'^Xr,s, 
which can of course be achieved by rescaling the couplings gr^s- In both cases the coupling g 
disappears from the expressions of the CF's. 

The Landau-Ginzburg formalism allows us to represent in a condensed form the essential 
properties of a topological field theory by means of a suitable (usually, polynomial) potential 
W. In the case of the topological field theory (let us call it this way) in question we must 
somehow expect an extreme situation. The reason is that Sq looks very much like an extension 
of the Ak type LG models to a model with an infinite number of fields in both positive and 
negative direction. So it seems natural to look for a LG model in which the ring of primary 
fields TZ is the ring of all monomials in both negative and positive powers of a certain variable ( 
and Q is represented by the identity. Let us call the latter ring V. Since we expect 7^ = V/dW 
we see that W can only be a constant at the critical point (in particular W will not be singular). 

We have such a potential at hand. Let us define 

W = LL (8.6) 

At the critical point we have W = x. Next we set 

dW dW 

</,,,, ^ ^^(^50) = M{r,s)x'C-^ </'o,o = ^(^o) = 1 

and in particular 
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We define the following map 

(l>n ^ Tn, 00 ^ Q, <P-n ^ T-n 

and claim that it is an isomorphism between the algebra TZ and the the field algebra A. It is 
elementary to prove the isomorphism by checking the few non-trivial cases. 

On the basis of the above argument we should not expect the formula for the three-point 
function in terms of the potential W to be exactly the same as in the ordinary LG models. In 
fact it is not, but it is very close to it 

< 4>ri,si(pr2,S2(Prs,ss >0= f "7 :^ (8.8) 



This reproduces exactly formula ( 3.19| ) with the obvious correspondences 4'ri,Si -^ Xri,Si, i = 
1,2,3. 
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